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ABSTRACT 

The adiabatic guiding center motion of an energetic charged particle in an 
idealized representation of the interplanetary magnetic field is discussed. Ex- 
pressions are presented for all drift velocity components. A comparison of 
these drifts establishes that, for the cases of physical interest, the electric drift 
velocity is by far the dominant drift velocity. 

Considering only the electric drift velocity, an expression is derived for 
the time of flight of a particle leaving the sun. Similarly, it is shown analytically 
that the particles remain attached to the same field line on which they left the 
sun; thus indicating their point of origin. High energy unscattered particles are 
shown to be rapidly channelled down the field lines, resulting in a highly aniso- 
tropic particle distribution in the vicinity of the earth. These analytic results 
are in excellent agreement with numerical integration of the complete drift 
equation and compare favorably with the experimental data. 
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ON THE DYNAMICS OF A CHARGED PARTICLE 
IN A SPIRAL FIELD 


I. INTRODUCTION 

Analytic solutions are presented which describe the adiabatic guiding center 
motion of an energetic charged particle in the idealized spiral model of the 
interplanetary magnetic field ( Parker , 1958). These solutions have simple 
physical interpretation and show the dependence of the motion on the various 
solar wind and particle parameters. Illustrative figures drawn for typical values 
of these parameters are also displayed. 

To estimate the validity of the approximate method used herein, we solved 
numerically the complete guiding center equation of motion. Our analytic solutions 
were confirmed in all cases by these computer results and agreed with most of 
the numerical findings of Winge and Coleman (1968). In addition, the analytic 
approach provides a generalization of the computer results. 

The solutions apply to all charged particles that are not appreciably scattered 
by magnetic irregularities in the interplanetary field. While many experimental 
results can be understood from the analysis herein, this paper is not meant to 
describe the exact motion of a charged particle from the sun to earth; disagree- 
ments between our results and experiments may arise because of the more com- 
plex and irregular structure of the actual magnetic field. 

H. MODEL OF THE SPIRAL FIELD 

We use the same model (Parker, 1958) for the interplanetary magnetic field 
as was discussed in some detail by Winge and Coleman (1968). The sun’s surface 
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magnetic field is combed out into interplanetary space by the radial expansion of 
the perfectly conducting solar wind. Since the magnetic flux lines lie along solar 
wind streamlines when viewed in a coordinate system rotating with the sun, the 
shape of the field is determined beyond r o , the radius at which corotation ceases. 

The magnetic field is shown in Figure 1 and is written in an inertial, helio- 
centric coordinate system as 

B = B o ( r </ r ) 2 [e r - fir sin e/W p ej (1) 

where the e 's are the unit vectors in the (r, 6 , 0) spherical polar coordinate 
directions with 6 the colatitude measured from the sun's rotation axis and 0 the 
longitude, which is positive in the direction of the sun's rotation. B o is the radial 
component of the field at r = r o , Q is the angular velocity of the sun and V p is 
the velocity of the solar wind, which is assumed constant and radial. For the re- 
mainder of this paper in numerical examples, we shall choose r o = 1.38 x 10 11 cm 
or twice the radius of the sun, B o = 5 x 10 -5 (lAU/r o ) 2 gauss = 0.591 gauss and 
fi= 2.90 x io“ 6 rad-sec -1 . V p shall be a constant parameter in the problem. 

The magnetic field may be simply expressed by using a unit vector e x , 
which lies along the magnetic field, as 

B = Bej = B o (r o /r ) 2 sec a e x (2) 

where 

e x = cos a e f - sin a (3) 

The angle a between B and e r is given by 
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a 


tan -1 (r/a) 


( 4 ) 


- tan" 1 (-B/B r ) = 

with the quantity a being a convenient scaling length: 


a - V p /(fi sin (9) (5) 

It is the radius at which the magnetic field makes an angle of 45° with respect to 
e r . For V p = 5 x 10 7 cm-sec' 1 , a = 1.72 x 10 13 cm at 6 = 90°. 


In the inertial system there is an electric field. 


E 


~V x B - ~(l/c) fi sin 9 B r 2 /r 

CP \ / oo & 


l 

-V p r B cos a/a e 0 


( 6 ) 


where c is the velocity of light. This electric field results from the convection 
of the magnetic field lines out with the solar wind speed. It should be noted that 
no electric field would be seen in a reference system moving with the solar wind 
plasma or in a reference system rotating with the sun. 


III. CHARGED PARTICLE DYNAMICS 

To simplify the discussion of the motion of a charged particle in a non- 
uniform electromagnetic field we use guiding center theory, that is, we concern 
ourselves with the motion of the guiding center, the point about which the particle 
may be considered to be gyrating. The particle itself may be shown to be nearby. 
The conditions under which guiding center theory is valid have been discussed by 
Alfvdn and Falthammar (1963) and presented for our problem by Winge and 
Coleman (1968) as the criteria for adiabatic invariance. These criteria were 
shown to be satisfied for all but extremely high energy particles. 
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The guiding center velocity V may be written in terms of a velocity V |( 
along the field and a velocity V d , the drift velocity, perpendicular to the field, as 

V = V„ + v d (7) 

The actual particle velocity is V plus V c , the velocity due to the gyration motion 
about the guiding center. 


Northrop (1963) has expressed the drift velocity of a particle with mass m, 
charge e and magnetic moment ^ as 


V 


d 



uc „ me 
cE + — VB + — 
e e 



dU E 
dt + dt 



where the electric drift velocity 


( 8 ) 


U E - c E x e j/B 
and the gravitational acceleration 

g = - M e G e r /r 2 

M 0 and G being the solar mass and the universal gravitational constant, respec- 
tively. The first two terms of the right hand member of Equation (8) are readily 
identified as the electric drift and the gradient B drift, while the remaining 
three terms account for the acceleration drift. 

Some care must be taken in the evaluation of the magnetic moment /j., which 
is an adiabatic invariant of the motion. The Appendix discusses this point and 
shows why for most cases we choose 
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( 9 ) 


m V. 2 sin 2 i p m V c 2 

^ “ 2B “ 2B 

where V T is the total particle velocity and 0 is the pitch angle, that is, the angle 
between the particle's total velocity vector and the magnetic field. 

IV. DRIFT VELOCITIES IN A SPIRAL FIELD 
A. Expressions for Velocities 

Using the representation of the spiral field given by Equations (1) and (6), 
we will now obtain an explicit formula for each drift term in Equation (8). This 
will permit us to demonstrate that the electric drift is by far the largest term, 
which will enable us to give a complete and accurate, although approximate, 
description of the guiding center motion. 


The electric drift is 


U E = cE x e x /B = 


V r 

p 


r 2 + a 2 


( r e r + 3 e rf>) 


( 10 ) 


The gradient of the field magnitude is 

VB = [(r 2 + 2 a 2 ) e 

r ( r 2 + a 2 ) L V 

which gives the gradient B drift as 


- .2 


r cot a e , 


( 11 ) 


U„ = -rr- x VB = (r 2 + a 2 ) -3/2 r 2 cot#e + (r 2 +2a 2 )e 0 + ra cot 6e± (12) 

“BVc/ g L r 


For a steady state problem, the operator in Equation (8) reduces to 
lowest order to the convective derivative along B ; 
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Hence the acceleration terms are 


de 


-Ai e B x V " dt eB 


“ = V„ 2 (r 2 + a 2 ) -372 (-r 2 cot 8 e r 


( 13 ) 


+ (r 2 + 2 a 2 ) e e - ra cot 8 


me e i dU E 


-A2 e B x dt eB 


me V p V u a 


(r 2 + a 2 ) 2 


(r 2 cot 8 e f - a 2 e g + ra cot 8 (14) 


Ug 


me M 0 G 


e B 


x 9 


e B r (r 2 + a 2 ) 172 


(15) 


We can combine results (10) through (15) by defining the total drift velocity 


as 


= V 


d r 


+ v 


ae e 8 


+ V 


d<£ e 4 > 


(16) 


where 


'dr 


V r 2 m c V 2 

P ^ r 2 (r 2 + a 2 ) 372 cot 6 — r 2 (r 2 + a 2 ) 372 cot 8 


r 2 + a 2 e 


m c v p V„ a r 2 
eB 


(17) 


(r 2 + a 2 ) 2 cot 8 


V = V dr a/r 


(18) 


V d0 = ^ ( f2 + 2a2 ) ( r2 + a 2 )"" 372 + 


m c V, 2 
e B 


(r 2 + 2a 2 ) (r 2 + a 2 ) 372 


(19) 


mcV V„ a 
eB 


(r 2 + a 2 )’ 2 - 


mcM. G 


e B r 


0 / 2 x 2 t/2 

(r z + a z ) 
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It is of interest to note that V dr and V d( ^ for particles in the equatorial 
plane are due entirely to the electric drift and that all other drifts are propor- 
tional to the mass of the particle. 

B. Comparison of Velocities 

To demonstrate that U E is the dominant drift, we proceed in the following 
manner: 1.) we find the largest contributor to the 6 drift; 2.) we show that the 
(r, 4>) drifts are at most comparable to their 6 counterparts; and 3.) we com- 
pare the term in part 1.) to U E and show it is small. 

1. ) We first consider the drift velocity in the 6 direction, given by Equa- 
tion (19). 

We compare the various terms to the second term U Aie : 

^ = - V aVfv,, (r 2 + 2a 2 ) (r 2 + a 2 ) 17 

J A10 L 

U Gg _ _ (r 2 + a 2 ) 

U Ai<9 r V,, 2 (r 2 + 2a 2 ) 

It may easily be shown from these expressions that U A10 is the dominant 6 drift 
except when the initial pitch angle is near 90° . 

2. ) We will now compare the various terms of the total (r, 0) drift 

■ [ v , 2 + v d 1/2 

with their corresponding 6 drifts: 


U 


b e 


V. 


U 


Md 


and 
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r cot 6 (r 2 + a 2 ) 1/2 /(r 2 + 2a 2 ) 


U b : 


U Ar V ,H e = - r cot 0 (r 2 + a 2 ) 1/2 /(r 2 + 2a 2 ) 


U A 2 " W t4> /V e ' ~ r cot 6 (r 2 + a 2 ) 1/2 /a 2 


U G • %/v, = 0 


These ratios are all of order lor less since r is of the same order as a . 


3.) Comparing the total electric drift velocity U E to U Al(9 , we find 


U Ai* _ mcV n /r 2 + 2a 2 ' 
U e = e B V p r r2 + a2 


( 20 ) 


To a good approximation V M is constant. Thus for V p = 5 x 10 7 cm-sec~ 1 
and protons with V n = 3 x 10 9 cm-sec _1 


U 


A10 




5 x 10" 20 


a r (2a 2 + r 2 ) 
(a 2 + r 2 ) 372 


This is a maximum at r = 2 M a and there is of order 10 -6 . For cases in 
which the initial pitch angle is close to 90°, U B(5 will be the largest contributor 
to 1.) when the particle is close to the sun. In this case V 2 would replace V 2 j 
in Equation (20) and the argument would proceed as before. 


Thus the electric drift dominates this problem and we are able to write that, 
to a very good approximation. 
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r V 


r 2 + a 2 


( re r + ae rf.) 


( 21 ) 


We note there is no 0 component. 

V. PARTICLE MOTION 

In this section we will derive expressions which describe the motion of a 
charged particle's guiding center in the spiral magnetic field. Specifically, we 
will consider an energetic particle (e.g., proton) of energy & and initial pitch 
angle «/» 0 which is injected from the sun into interplanetary space at r = r 0 , 
colatitude & 0 and longitude 4> 0 . Our task is to obtain r(t) and <£(t) as functions 
of these initial parameters and the solar wind velocity. Since V e % 0, 6 (t) = 8 0 , 
within our approximation. 

A. Radial Motion 

The radial position of the guiding center is the solution of the differential 
equation 

d r 


dt = V r = V dr + \ 


( 22 ) 


where 


V dr = V p r 2 / (r 2 + a 2 ) 


(23) 


and 


- V„ cos a - V n a/(r 2 + a 2 ) 1/2 


(24) 


We use conservation of energy and the adiabatic invariance of the magnetic 
moment to determine V n (r). The particle's kinetic energy is a conserved quan- 
tity because changes in the gravitational energy are entirely negligible and 
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because the electric potential energy is automatically a constant of the motion 
since the particle moves on a surface of constant 8 . If V T = (2£/m) H is the 
initial speed of the particle, then 

V T = V n + V c + V d ( 25 ) 

r\ 

V d 2 (r) is given by Equation (21) and V c 2 (r ) = B(r) from Equation (9). 
Using Equations (2) and (9) 


V c 2 (r) 




Thus Equation (25) reads 


(26) 





sin 2 x[r 0 



(27) 


Using this expression for V |( (r) in Equation (22) gives, in principle, r (t). Un- 
fortunately, we have not been able to perform the required integral, but some 
well founded approximations will enable us to arrive at an explicit representation 
of t (r). 

We first neglect V d 2 compared to V T 2 since V d < V p 2 < < for all particles 
under consideration. For example, a 1 Mev proton has V T = 1.38 x 10 9 cm-sec” 
while a typical solar wind velocity is 5.00 x 10 7 cm-sec -1 . Next, noting that 
r Q 2 < < a 2 for all 9, we can replace Equation (26) by the approximate form 

V 2 £ r 0 2 sin 2 i/^/r 2 (28) 
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Differences between this and the more exact Equation (26) are significant only 
for r > a , when V c 2 < < V. 2 . Thus 


V„ ~ V T [l - r 0 2 sin 2 i/< 0 /r 2 ] 1 


/ 2 


(29) 


Using similar arguments, we neglect V dr compared to V (| and obtain the 
radial equation of motion 


dr 

'v, V 

dt ~ v " r 


a 

= V Tr 


(r 2 - r 2 sin 2 Q 1/2 
(r 2 + a 2 ) 1/2 


(30) 


which may be integrated to give 


t ( r ) 


2a\~ |( r2 " r o 2 sin2 to) 1/2 O 2 + a2 > 1/2 " r o c °s <£o (r 0 2 + a 2 ) 


(31) 


+ (a 2 + r Q 2 sin 2 In 


(r 2 + a 2 ) 1/2 + (r 2 - r 2 sin 2 1 
( r o + a2 ) 1/2 + r o cos to 


This gives the flight time for a particle to reach a distance r from the sun if it 
is injected at r 0 with velocity V T and pitch angle i// 0 . It is plotted for various 
initial conditions in Figure 2. 

If we consider r Q < < r, a, then Equation (31) for 6 = 90° simplifies to 
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It is easily shown that for large r , t ~ r 3/2 , which agrees with the asymptotic 
slope of 3/2 seen in Figure 2. 

Lin and Anderson (1967) have tabulated over 25 prompt solar electron 
(>40 Kev) events for the years 1964-1966. The shortest travel time recorded 
for the electrons was 1500 seconds, while the longest was 5000 seconds. A 
travel time of 2000 seconds was typical. We find that 40 Kev electrons would 
move from r 0 to the earth in our model in 1950 seconds, 1550 seconds and 
1400 seconds for V p = 200 km-sec -1 , 350 km-sec -1 and 500 km-sec -1 , respec- 
tively. These results are nearly independent of the initial pitch angle and, since 
they deal with the lowest energy particle measured, they are the longest travel 
times that would be measured according to our theory. The fact that larger 
travel times are measured indicates that path length is longer than assumed (a 
more irregular field) and/or injection of the particles at the connecting field 
line takes place later in the event. The latter explanation would mean that the 
particles either take some time to diffuse to the correct field lines at the sun or, 
more simply, particle injection occurs after electromagnetic emission. 

Fan, Pick, Pyle, Simpson and Smith (1968) have reported solar proton (0.6 - 
13.0 Mev) events during 1966 and have found travel times of typically 10 4 
(± 2 x 10 3 ) seconds. We illustrate in Figure 3 how the time of flight of a 13.0 Mev 
proton in the simple spiral field would vary with the solar wind velocity V p from 
Equation (31). The flattening of the curve at high V p illustrates the straightening 
of spiral field lines. Once again, we emphasize that these results are a lower 
bound on the problem. The remarks of the above paragraph also apply here. 
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B. Azimuthal Motion 


A calculation of the azimuthal motion is most easily done in a reference 
frame which is rotating with the sun. If and denote the velocity components 
of the guiding center in this frame, then 


V M 3 V 2 1 

V' = V = + — = 

(r 2 + a 2 ) 1/2 (r 2 + a 2 ) ( r 2 + a 2 ) 1/2 


V„ a + 


V r 2 

p 


(r 2 + a 2 ) 172 


(33) 


v; = V r flr Sin 0 


V " r V p 3 r „ 

+ -fir sin 6 

(r 2 + a 2 ) 172 r 2 + a 2 


a(r 2 + a 2 ) 172 


V„ a + 


V* 


(r^ + a z ) 


(34) 


Clearly, 


v ; 

= - r/a = B^/B r (35) 


which shows that the guiding center trajectory moves out a spiral magnetic field 
line when viewed in the rotating (solar) frame. This occurs because there is no 
electric field in the rotating frame and, hence, no electric drift. An inertial 
observer would describe the motion as analogous to that of a bead sliding outward 
along a spiral wire, which is rotating with the angular velocity ft, or would say 
that the particle's <j> is the 4> of the spiral field line minus the <f> through which 
the sun has rotated during the particle's flight. This result may be expressed 
analytically as 

ft sin 8 ^ 

4> ~ 4> 0 ~ v [ r o " r ] + 0 t(r) (36) 

P 
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with t (r) being given by Equation (31) and 0 O being 0 at r 0 ; it is shown in 
Figure 4. We note that a solar wind particle, which is injected radially at r Q , 
with V T = V , has 0 - 0 O . 

The last term in Equation (36) is important if the particle energy is small 
or if the field is tightly coiled, that is, if V p is small or r is large. For most 
cases, however, we are justified in ignoring the last term as well as r 0 . Then, 
if a particle is received at earth, we can say, within our approximate model, 
that the flare causing the particle injection occurred at 

0 O = fir /V p West (37) 

This is an equation that has been used in the past but now it has somewhat more 
justification. Lin and Anderson (1967) have shown that the average heliographic 
longitude for simple (using their terminology) electron events is 58° West. This, 
from Equation (37) , would correspond to an average solar wind velocity of about 
425 km-sec -1 , quite a reasonable value for solar minimum. 

C. Particle Anisotropy 

We have shown in the previous section that particles move out the magnetic 
field lines upon which they were injected. As a last point, we show how the parti- 
cles are rapidly channelled down the magnetic field lines. The average pitch 
angle 0 seen by a stationary observer is given by 

sin 2 0 = (V 2 + V d 2 )/V T 2 (38) 

Using Equations (21) and (26), we may write 
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(39) 




sin 


-l 




sin 2 </r Q 


r 2 v p 2 /v T 2 

r 2 + a 2 



We see that at r = r 0 , <// % </> 0 and at r > > a, </» % sin" 1 (V p /V T ). The first term 
is energy independent and measures the effect of B(r) on V c (r). The complete 
behavior of ^ as a function of r for V p = 500 km-sec" 1 is shown in Figure 5, 
where the asymptote of sin -1 (V /V T ) is clearly small. Hence at the earth, 
any unscattered particles should have a highly anisotropic distribution; the 
anisotropy axis should lie near the field direction. These conclusions are in 
substantial agreement with the experimental results of Fan et al . (1968) and 
McCracken, Rao and Bukata (1967). 

VI. SUMMARY 

We have investigated the dynamics of an energetic charged particle in a 
spiral model of the interplanetary magnetic field. After showing that the electric 
drift velocity is the only large drift velocity, we were able to arrive at approxi- 
mate analytic solutions for the adiabatic motion of the guiding center of the 
particle. 


These solutions show that the guiding center of the particle moves like a 
bead sliding outward along a spiral wire which is rotating with the angular 
velocity of the sun. The spiral wire represents the field line on which the 
particle starts. A reasonably simple expression is presented for the time of 
flight, which is inversely proportional to the total particle velocity but which has 
a more complex dependence upon the distance from the sun, the solar wind 
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velocity and the initial pitch angle. Our calculations indicate that the particles 
have a highly anisotropic distribution upon arrival at the vicinity of the earth. 

The above findings are shown to be in substantial agreement with experi- 
mental data on unscattered solar cosmic rays. The expressions presented 
generalize and lend physical interpretation to the computer results of Winge 
and Coleman (1968). 
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APPENDIX : EVALUATION OF fj, 


The magnetic moment /a is a measure of the flux enclosed when the particle 
circles the guiding center. As an adiabatic invariant of the motion, it is a 
powerful tool in solving for the motion along the field lines. The classical 
description is 


fi - mV c 2 /2B - mV T 2 sin 2 \p/ 2B (A-l) 

The last step is only true when V d < < V T sin <//. At other times care must be 
exercised. We define y in terms of the instantaneous vector velocity perpen- 
dicular to the field line V T 

fi = "(V^- V d ) 2 /2B (A-2) 

The case of zero magnetic moment arises when the particle is injected with 
just the drift velocity (V I( = V d ). This is of particular interest; the particle 
does not gyrate at all, is unaffected by a gradient B drift and continues to move 
across the field at V d . This illustrates how a solar wind particle is permitted 
to move across the interplanetary magnetic field. 

There is some ambiguity in Equation (A-2) because Vj ^ is usually defined 
in terms of only an angle 0 and an energy. Thus the direction of V x ^ is un- 
known; we arbitrarily choose it to lie in the positive 4> direction in the (r, <p) 
plane. 

The magnitude of Vj ^ is V T sin i/>, which will dominate V d in this problem 
for all pitch angles except when ~ 0 (V d /V T ). It must be emphasized, however, 
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that this is a weak restriction. Thus for almost all cases Equation (A-l) is an 
accurate expression for the magnetic moment. 
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FIGURE HEADINGS 


Figure 1. Spiral field configuration, showing parameters and coordinate system. 

Figure 2. Time of flight versus r for protons injected at r„. 

Figure 3. Time of flight for a proton injected at r 0 to reach 1 A.U. versus the 
solar wind velocity V p . 

Figure 4. Actual proton trajectory. Dots show the position of the particle in 

4 

intervals of 10 seconds. 

Figure 5. Average pitch angle versus r for protons injected at r 0 . Note that i p 
decreases rapidly with increasing radius. 
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